General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 




( M AS A-CE- 174073) TUE flAfillAA CLIMATE: NeS-12e49 

EMEBGi BALANCE MODELS HITE CC2/H20 
ATHOSPHEBES Seaiaaaual Status fieport, 1 

Mar. - 1 Sep. 1984 (Netr ycxk Oniv., Nev Unclas 

Xork.) 11 p EC A02/ME A01 CSCL 03B G3/91 19781 











; . 


J V 




< ■■ V . . ■ . . . ^ • 







SEMI-ANNUAL STATUS REPORT 


THE MAFITIAN CLIMATE: ENERGY BALANCE MODELS 
WITH CO 2 /H 2 O ATMOSPHERr«J 

Principal Investigator 
Martin I. Hoffert 
Prepared for 
NASA GRANT NAGW-573 
March 1, 1984 - September 1, 1984 



NEW YORK UNIVERSITY 

FACULTY OF ARTS AND SCIENCE 

DEPARTMENT OF APPLIED SCIENCE 

26-36 Stuyvesant Street 
New YorK, N.Y. 10003 



The Martian Climate: Energy Balance Models with Carbon Dioxide/Water 
Atmospheres 

Progress Thus Far 

As described in our original proposal and work plan (including some 
modifications per the suggestions of the MDAP Review Panel), we are 
developing a multi-reservoir, time dependent energy balance climate 
model for Mars driven by prescribed insolation at the top of the 
atmosphere. The first approximately half-year of the program was devoted 
to assembling and testing components of the full model. Specific 
accomplishments have been made on a longwave radiation code, coupling 
seasonal solar input to a ground temperature simulation, and conceptualizing 
an approach to modeling the seasonal pressure waves that develop 
in the Martian atmosphere as a result of sublimation and condensation of 
COg in polar regions. These steps are described below. 

(1) Longwave radiation code 

The previous study by Hoffert et al. (1981) of Mars' surface 
temperature distribution as a function of carbon dioxide and 
water vapor atmospheres utilized a 3-band longwave radiation code 
for the outgoing IR flux from the planet, i.e., the top-of-the- 
atmosphere (TOA) flux. The seasonal energy balance model (EBM) 
under development with separate atmosphere and surface temperatures 
requires the backradiation from the atmosphere to the surface at 
the bottom of the atmosphere, i.e., the BOA flux. A number of 
model runs using the longwave code integrated toward the surface 
have been made to develop BOA correlations as a function of 
surface pressure and temperature for use in the mul n -reservoir EBM. 
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(2) Seasonal Ground Temperature Simulations: Current Mars 


Expressions for the seasonal and latitudinal variation of 
solar radiation at the top of the atmosphere have been 
formulated as a function of orbital radius, eccentricity, 
obliquity and longitude of perihelion. This will allow 
consideration of the effects of astronomical long-period 
variations in orbital elements on insolation, and hence, 
on climate. This has been used to drive a "forcu-restore" 
model for ground temperature (neglecting at this point the 
influence of atmospheric and boundary layer heat transfer). 

The method is derived in Hoffert and Storch (1979), is 
computationally faster than a point-by-point integration 
of the ID heat conduction equation into the ground; it 
involves a numerical solution of a ordinary (as opposed to 
a partial) differential equation: 

(l/2)C5dT^/dt= s(^()i)(l-a)- coT* - ( 1 /2) Cj!2( T 5 -<Tj> ) , (1) 
where 

Cg=(2/n)I, where I is the thermal inertia 

Tg, <Tg> are the ground surface temperatures, instantaneously 

and annually averaged, respectively 

t = time 

a = solar albedo 
e = longwave emissivity 
= angular velocity of periodic forcing . 

Using e = 1.0 and I = 272 J/m^-K-s^^^, the seasonal variation of 
daily mean surface temperatures at latitudes of the VL-1 
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and VL-2 lander sites are shown and compared to lander 
anemeometer temperature data (from Seiff, 1982) in Figs. 

1 and 2. The thermal inertia is within the range of Mars 
values I 'V' 70-200 J/m^-K-s^^^ cited by Jakosky (1979). These 
Mars surface values of I were derived by curve fits to a full 
integration of the unsteady heat conduction equation (in depth 
and time) by Kieffer et al. (1977). Equation (1) is exact for 
T^ for harmonic (sinusoidal )forcing and approximate when the 
forcing is periodic but nonharmonic; it avoids the numerical 
integration of PDEs at each surface point (an advantage for 
global climate models) and is more accurate than a slab model, 
since it calculates the temperature at the surface, as opposed 
to a mean slab temperature. 

Fig. 1 shows the model versus data for the VL-1 site at 48°N 
fora=0.25. The agreement is generally quite good, with a 
seasonal amplitude of about 60°K. The model precedes the data 
in phase by about 10-20° of aerocentric longitude, indicating 
the possibility that I may be greater for this site than the 
Martian average. This can undoubtedly be improved by tuning 
the local I as described below. 

Fig. 2 compares the model against data for the VL-1 site at 
22°N, using a =0.32. Here the agreement is again quite good 
except in the vicinity of a temperature minimum during northern 
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hemisphere winter near 300° aerocentric longitude. This is 
likely due to an insolation cutoff by the dust storm occuring 
between 275°< Lg<320°, which "thickened" the atmosphere to an 
optical depth of 8 (James and North, 1982). Notice that 
= 275° is where the model VL-1 temperature falls' below the 
model. There v/as also a dust storm between 205°< Lg<250° , 
but its effects on atmospheric transmissivity was about two - 
orders of magnitude less than the storm at the temperature 
minimum. The dust storm effects do not show up in the model 
versus data comparison of Fig. 1; either the dust loading was 
less further north, or since the temperature minimum of about 
157°K was near the condensation frost point or C02» latent 
heat effects prevented the temperature from falling any lower. 

This is also to be investigated further. 

To exploit Kieffer's (1977) computations of I and a for Mars, 
we have obtained from H. Kieffer the Mars Average Data Set (MADS) 
containing temperature and albedo data. Our intention is to check 
the seasonal surface temperature model against data for more 
latitudes than the two Viking lander sites. However, there are 
significant amounts of missing temperature data (Martin, 1981), 
and it is not yet clear the extent to which valid daily T^ values 
can be obtained for all latitudes and seasons. This investigation 
is underway. Nevertheless, it appears we can utilize the albedo 
data from MADS (for non-dusty times) as a function of latitude and 
possibly time of year. as an input to our seasonal model, or 
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possibly parameterize seasonal albedo change as a function of 
other properties (such as T^). The distribution of I as a function 
of latitude is being obtained from the full Mars Consortium Data Set. 

(3) Seasonal Pressure Waves 

Despite the reasonable (present day) surface temperatures at 
middle latitudes neglecting atmospheric transport, it is necessary 
to include meridional mass transport of CO 2 in our seasonal 
model to accoiTiC for seasonal pressure waves associated with 
CO 2 cap sublimation/condensation cycles. The James and North 
(1982) seasonal model assumes PCO 2 globally equilibrates over 
the duration of their numerical timestep of 'v 4 sols. Phil 
James (personal communication) cittis Haber le's Mars GCM results 
indicating a relaxation time ( by sublimation-driven winds) from- 
the north pole to VL-1 as no more than 'V/ 8 sols. • However, the 
pressure data for VL-1 and VL-2 (Sieff, 1982) show a minimum 
pCO^ at VL-1 (22°N) during northern hemisphere winter, when 
CO 2 is condensing into the southern polar cap, to precede the 
corresponding minimum at VL-2 (48°N) by perhaps 0-30 sols. 

Similarly, the peak at VL-1 during southern hemisphere summer, 
when mass flux of CO 2 is entering the atmosphere from sublimation 
at the south polar region, precedes the maximum at VL-2.. ■ 

To study this further, we are obtaining VL-1 and VL-2 pressure 
data from J. Tillman (2245 sols @ VL-1 and 1050 sols 0 VL-2-- 
He expects to release this data to NSSDC in December, and we 
should receive it shortly). Our intent is to use it to estimate 
the meridional seasonal mass flux of CO 2 by examining the existence 
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of phase lags between VL-1 and VL-2, and using their presence 
to calibrate equator- to-pole mass flux rates in the seasonal 
model. These considerations apply to meridional fluxes within 
the northern hemisphere only. Another effect which may impact 
the observed seasonal pressure waves is blocking by the 
intertropical convergence zone (ITCZ). Both Haberle's model 
and Kahn's data show a well -developed Hadley circulation in 
the tropics which could well effect interhemispheric transports. 
These effects, individually or in combination, could well effect 
the observed seasonal pressure waves, and will be studied in 
relation to our modeling of the current Martian climate. 


- 6 - 




TEyPCRATURE 00 


ORlQiNAL PAGtl 
OF POOR QUALITY 



aerocentric longitude (DEGREES) 


figure 2 


R»f*r*ncM 


Koff* *t» n.l.a and J. Btoreli (i?7?> A sch«a* for computinf^ sur/aeo 
ClusM from noon tlo«# obsorvations . » BOUNDARY LAYER METBOROLO IY» vol 
17» pp.429— ■ D.- Raidal Publishing* Dordrocht* Holland 

I 

Hoffort* M.I.* A.J. Callogari, C.T. Haloh* and V. Zloglor <1981) 

Liquid wator on Mars: an snsrgy balancs modsl for C02/H20 
af '^^sphsros. ICARUS* vol. 47* pp. 112-129 

Jakosky* B.M. (1979> Ths affects of nonideal surfaces on the derived 
thermal properties of Mars.* J. CEOPHYS. RES.* vol. 84* pp. 

82S2-8262. 


James* P.B.* and C.R. North (1982) The seasonal C02 cycle on Mars: an 

T 

application of an energy balance climate model* J. CEOPHYS. RES.* 
vol. 87* pp. 10271-10283. 

t 

Martin* T.Z. (19811 Mean thermal and albedo behavior of the Mars 

• ■ 

surface and atmosphere over a martian year.* ICARUS* vol. 45* pp. 
427-448 

Seiff, A. (19B2> Post-V*king models for the structure of the summer 
atmosphere of Mars.* ADV. SPACE. RES.* vol. 2* pp. 3-17. 


Ui\iw, 


Q£ £OOR QUAiunX 


OKJUR7AE PAXTE IS 

OF POOR 



